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INTRODUCTION 


This  report  covers  special  work  carried  out  at  the  University  of  Arizona, 
Optical  Saences  Center  under  support  provided  through  the  Joint  Services  Optical 
Program  (JSOP). 

The  work  summarized  in  the  report,  and  detailed  in  the  appendices,  is 
responsive  to  requests  from  two  sponsors.  In  each  case  the  work  required  system 
analysis  and  technical  evaluation  of  the  potential  in  two  areas  of  optical 
technology.  The  work  on  a  high-speed  video  system  was  led  by  Professor 
Eustace  L.  Dereniak.  The  work  on  membrane  mirros  was  led  by  Professor 
William  L.  Wolfe. 


Robert  R.  Shannon 
Director 
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HIGH-SPEED,  HIGH-RESOLUTION  VIDEO  IMAGING 
SYSTEM 


E.L.  Dereniak 


SCIENTIFIC  PERSONNEL 

E.L.  Dereniak 
T.  Graeve 
L.  Flath 


OBJECTIVES 

Our  objectives  were  to  study,  design  and  implement  test  instrumentation 
for  high-speed  imaging  devices  and  instrumentation  for  obtaining 
photometric/radiometric  data  needed  to  evaluate  the  in-flight  performance  of 
such  devices. 


RESEARCH  FINDINGS 

1000s  Imaging  System  Study 

This  part  of  the  project  involved  the  study  of  strategy  and  specifications  for 
a  ruggedized  1 000-frame-per-second,  airborne  qualified,  megapixel  (1024-by-1024) 
imager  and  support  system.  The  initial  effort  of  our  research  on  high  speed 
visible  solid  state  detectors  was  to  theoretically  address  the  radiometric  aspects 
of  the  problem  to  estimate  the  performance  of  a  particular  technology 
configuration  with  varying  figures  of  merit 

In  our  initial  investigation  we  used  a  radiometer,  a  filter  and  some  targets 
to  substantiate  our  theory  with  experimental  data.  The  supporting  calculations 
and  the  results  are  shown  in  Appendix  A.  The  solar  constant,  as  well  as 
radiometric- to-photometric  conversion  factors  were  calculated  for  various  spectral 
bands.  The  filter  used  in  the  experiment  transmitted  from  036  \an  to  0.66  jim. 
The  projected  signal  values  for  this  spectral  band,  given  a  standard  CCD  imager 
of  (10  jimP  pixel  size,  25%  quantum  efficiency  and  1  ms  integration  time,  range 
from  10?  to  10^  electrons/pixel  for  scenes  illuminated  by  direct  sunlight 
Appendix  A  charts  measured  vs.  calculated  irradiance,  using  standard  reflectance 
values  for  some  of  the  measured  targets.  The  close  agreement  between  the 
measured  and  expected  values  confirms  our  model. 
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Next,  we  analyzed  the  signal-to-noise  characteristics  of  a  standard  CCD 
imager.  The  data  and  calculations  are  shown  in  Appendix  B.  As  is  evident  from 
the  plots;  signal-to-noise  ratios  of  10  and  greater  can  be  easily  achieved  for  a 
variety  of  signal  and  read  noise  values. 

Appendix  C  shows  the  results  of  the  power  dissipation  study  for  a  frame 
transfer  CCD.  The  array  modeled  is  the  Ford  Aerospace  F1024B  10242  pixel 
imager.  Most  of  the  power  in  the  array  is  generated  by  capacitive  charging  and 
discharging  of  the  clock  lines.  Using  our  model,  the  total  power  dissipated  by 
this  array  running  at  1000  frames  per  second  was  estimated  at  approximately  10.5 
watts.  This  shows  the  need  for  reducing  the  line-to-line  capacitance  in  the 
polysilicon  clock  lines,  and  maybe  even  strapping  the  clock  lines  using  aluminum 
strips  to  bring  the  clock  signal  out  onto  the  device.  The  dissipated  power  can  be 
reduced  to  approximately  1  watt  just  by  taking  simple  measures  to  decrease  the 
line-to-line  capacitance. 

We  performed  additional  modeling  on  the  effects  of  the  frame-transfer 
process  on  image  contrast.  During  readout,  each  pixel  collects  additional  charges 
as  it  is  moved  across  the  device.  These  charges  reduce  the  contrast  in  the  image 
by  averaging  or  smoothing  the  pixel  contents  in  the  readout  direction.  For  a  full 
frame  architecture  device,  this  can  lead  to  a  contrast  reduction  of  50%.  In  a  frame 
transfer  device,  a  contrast  of  around  90%  can  be  expected.  In  an  interline  sensor, 
this  effect  is  negligible.  The  calculations  and  results  are  given  in  Appendix  D. 

Finally,  we  modeled  the  Ford  F1024B  device  on  SPICE  to  determine  the 
clock  pulse  degradation  due  to  the  polysilicon  resistance  for  high  frame  rates. 
The  model  showed  significant  exponential  rise  and  decay  effects  on  the  pulse 
shape  at  the  center  of  an  array  clock  line.  These  effects  are  due  to  the 
combination  of  the  poly  line  resistance  and  the  large  capacitances  between  the 
clock  line,  adjacent  lines  and  the  substrate.  The  maximum  working  array  clock 
speed  for  the  F1024B  device  is  500  kHz.  The  plots  shown  in  Appendix  E  further 
explain  the  model  used  and  the  resulting  pulse  shape  degradation. 


Airborne  Multispectral  Radiometer 

As  part  of  the  effort  to  further  analyze  the  working  environment  of  a  10003 
imaging  system,  we  designed,  built  and  tested  an  airborne  multispectral 
radiometer.  The  primary  function  of  this  instrument  is  to  measure  radiation 
reflected  off  an  airplane  in  several  different  spectral  bands.  These  data  then  will 
be  used  to  develop  exact  specifications  for  a  high-speed,  high-resolution  imaging 
system  recording  a  similar  scene. 
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The  system  is  based  around  the  Xybion  MSC-02  Multispectral  Video 
Camera  manufactured  by  Xybion  Electronic  Systems  of  San  Diego,  California. 
The  camera  uses  a  Sony  interline  transfer  CCD  with  786  (H)-by-493  (V)  pixels  and 
8.8  mm-by-6.6  mm  imaging  area.  The  CCD  is  located  behind  a  rotating  filter 
wheel  equipped  with  six  narrow  bandpass  interference  filters  centered  at  400  nm, 
450  nm,  500  nm,  600  nm,  700  nm  and  800  nm.  The  filter  wheel  rotates  at  600  rpm 
and  is  synchronized  with  the  RS-170  video  output,  so  that  each  consecutive  video 
field  corresponds  to  a  different  spectral  band. 

The  Xybion  camera  has  been  modified  to  include  two  photodiodes 
mounted  behind  the  filter  wheel  at  120°  to  the  CCD.  The  field  of  view  of  one  of 
the  photodiodes  has  been  matched  to  that  of  the  camera,  whereas  the  other 
photodiode  has  an  extremely  narrow  field  of  view  for  spot  radiometry.  The 
signals  from  the  photodiodes  are  processed  and  converted  into  digital  format 
using  two  pairs  of  logarithmic  amplifiers  and  16-bit  A/D  converters,  and  then 
time-delayed  for  a  fixed  number  of  fields  to  account  for  the  spectral  offset 
between  the  photodiode  data  and  the  camera  video  output.  The  data  then  are 
ready  to  be  sent  out  via  the  parallel  interface  to  the  Video  Data  Encoder 
(manufactured  by  Datum  Inc.  of  Anaheim,  California)  which  inserts  it  into  the 
current  video  field. 

The  multispectral  radiometer/camera  is  mounted  inside  a  heavy-duty 
aluminum  chassis  that  is  designed  to  withstand  the  harsh  conditions  of  external 
mounting  on  an  aircraft.  This  chassis  was  designed  at  the  University  of  Arizona 
and  was  built  at  Eglin  Air  Force  Base  under  the  supervision  of  D.  Snyder  and  H. 
McCormick.  It  includes  a  power  supply,  custom  electronics,  the  Xybion  camera, 
and  the  optics  for  the  radiometers  and  the  camera.  A  limited  amount  of  space  is 
available  for  additional  cameras.  The  unit  is  sealed  against  the  environment  with 
cork  gaskets,  but  there  are  valves  mounted  in  the  front  cover  plate  to  equalize 
pressure  differentials  between  the  inside  of  the  box  and  the  outside  air  as  the 
aircraft  changes  altitudes.  All  electrical  connectors  are  easily  accessible  on  top  in 
the  rear  cover  plate.  Appendix  F  shows  a  sketch  of  the  multispectral 
radiometer/ camera,  as  well  as  a  diagram  of  the  filter  wheel  and  a  list  of  the 
system  specifications. 


High-Speed,  High-Resolution  CCD  Test  Station 

To  evaluate  CCD  imagers  developed  for  use  in  the  10003  imaging  system 
described  above,  a  high-speed  focal  plane  array  test  station  was  built.  This  test 
station  consists  of  a  HP  Vectra  486  computer  for  waveform  generation,  data 
acquisition  and  data  analysis/display,  as  well  as  custom-built  hardware  for 
addressing  the  FPA  and  reading  out  the  data.  The  system  is  capable  of  running 
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multi-channel  arrays  (up  to  64  output  channels)  at  pixel  rates  up  to  one  MHz. 
Figure  1  shows  a  block  diagram  of  the  facility. 

The  custom  hardware  necessary  to  operate  the  test  station  consists  of  clock 
drivers,  bias  regulators,  pre-amps,  a  64-channel  multiplexer  and  a  correlated 
double  sampling  (CDS)  unit.  It  has  been  optimized  for  running  the  Ford 
Aerospace  F1024B  CCD.  The  clock  driver  and  bias  regulator  hardware  is  capable 
of  separately  adjusting  the  active  area,  storage  area  and  serial  register  clock  levels 
between  -10  volts  and  +10  volts.  The  DC  bias  levels  can  be  adjusted  between  -35 
volts  and  +35  volts.  For  output  signal  processing,  12  boards  containing  64  pre¬ 
amps  have  been  designed,  assembled  and  tested.  The  64-channel  multiplexer  card 
with  built-in  gain  and  offset  controls  is  used  to  switch  the  64  video  signals  into 
one  single  channel.  Finally,  the  CDS  unit  processes  the  video  signal  and  reduces 
noise  introduced  by  capacitive  elements  in  the  circuit.  It  also  includes  the 
necessary  control  circuits  for  fine-tuning  the  timing  pulses  for  sampling  and  A/D 
conversion. 

The  test  head  for  the  focal  plane  array  consists  of  a  four-layer  PCB  holding 
a  pin  grid  array  socket  into  which  the  117-pin  chip  carrier  for  the  CCD  is  inserted. 
The  PCB  is  mounted  inside  a  sealed  chamber  on  top  of  a  5-watt  thermoelectric 
cooler.  The  chamber  can  be  evacuated  to  eliminate  condensation  of  moisture  on 
the  cooled  device.  Electrical  connections  for  inputs  and  outputs  to  the  device  are 
made  through  three  32-pin  sealed  Bendix  connectors.  The  test  chamber  cover 
includes  a  2-in  window  and  an  adjustable-focus  C-mount  lens  holder.  Three 
thermistors  are  mounted  inside  the  test  chamber  for  temperature  monitoring. 

Several  pieces  of  software  have  been  developed  for  the  test  station.  The 
clock  waveforms  are  generated  by  a  data  generator  board  plugged  into  the 
computer  bus.  It  has  been  programmed  to  run  the  test  array  in  frame-transfer 
mode  at  speeds  ranging  from  <1  frame  per  second  (fjps)  up  to  100  fps.  Most  of  the 
testing  has  been  performed  at  either  6  fps  or  30  fps. 

A  1  MHz  ADC  board  is  used  to  acquire  the  data  presented  by  the  CDS 
unit.  Software  has  been  written  both  to  acquire  and  display  single-channel  data 
and  64-channel  composite  imagery.  In  single-channel  mode,  the  test  station 
retrieves  data  only  from  one  of  the  64  channels  of  the  device  and  displays  it  on 
the  screen.  In  64-channel  mode,  the  acquisition  hardware  automatically  switches 
to  the  next  channel  after  every  frame,  thus  assembling  a  complete  512-by-1024 
pixel  image  in  64  frames.  Both  modes  feature  a  real-time  display  with  256  grey 
levels. 
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temperature  readout 


Figure  1.  Block  diagram  of  100  fps  test  station. 
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RESULTS  AND  RECOMMENDATIONS 

Test  procedures  to  conduct  mean-variance  analysis,  quantum  efficiency, 
charge-transfer  efficiency  and  MTF  tests  have  been  studied.  Both  software  and 
hardware  to  conduct  automated  mean-variance  testing  has  been  developed.  The 
results  for  the  mean-variance  tests  thus  far  have  been  limited  by  a  very  high  read 
noise  level  in  the  acquired  data  (on  the  order  of  1000  electrons).  This  is  expected 
to  decrease  as  the  electronics  are  further  refined.  The  mean-variance  analysis  is 
further  complicated  by  the  fact  that  the  conversion  gain  in  the  system  is  kept  low 
to  stay  within  the  voltage  swing  of  the  electronics.  This  leads  to  a  low  conversion 
gain,  which  makes  the  mean-variance  curve  difficult  to  interpret. 

The  large  voltage  swing  in  the  output  signal  is  caused  by  a  DC  offset  of 
several  volts  between  the  reset  and  the  signal  levels  in  the  waveform  coming  out 
of  the  device.  This  DC  offset  is  caused  by  capacitive  effects  as  the  floating 
diffusion  on  the  chip  is  released  from  the  reset  potential.  The  DC  offset  falls  on 
top  of  the  signal  since  the  reset  clock  also  transfers  the  signal  charge  into  the 
floating  diffusion.  An  additional  output  transfer  clock  on  the  device  would 
eliminate  this  problem  by  providing  an  intermediate  sampling  level  to  be  used 
as  a  reference  in  the  correlated  double  sampling  circuit  The  electronics  could 
then  be  timed  to  amplify  the  actual  signal  and  ignore  the  DC  offset. 

Another  aspect  of  the  device  that  has  been  studied  is  the  variation  in  signal 
between  different  outputs  on  the  array.  Both  shifts  in  the  DC  offset  described 
above  and  in  output  sensitivity  have  been  observed.  The  output  sensitivity  drops 
by  up  to  three  orders  of  magnitude  towards  the  center  of  the  device.  Outputs  at 
the  edge  of  the  device  are  most  sensitive  and  those  at  the  center  of  the  array  have 
the  least  sensitivity.  Further  study  is  necessary  to  find  the  exact  cause  of  this 
effect,  but  it  is  believed  to  be  related  to  a  degradation  in  clock  or  bias  voltages 
towards  the  center  regions  of  the  CCD.  To  counteract  the  signal  variations,  a 
two-point  software  correction  scheme  has  been  implemented  to  equalize 
differences  in  both  offset  and  gain  between  the  channels.  However,  the 
magnitude  of  the  variations  exceeds  the  range  of  the  A/D  converter  used  for  the 
data  acquisition.  The  correction  should  be  performed  in  hardware  on  the  analog 
signal  to  work  satisfactorily. 
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W.L.  Wolfe 

SCIENTIFIC  PERSONNEL 

W.L.  Wolfe 
L.  Brod 
D.  Eklund 
V.  Sinclair 

OBJECTIVES 

Our  objectives  were  to  study  th in-film  mirrors,  under  vacuum,  whose 
properties  have  been  directionally  altered  by  x-rays  or  other  means. 


RESEARCH  FINDINGS 

Our  collaborator.  Dr.  William  Gamble  of  the  U.S.  Army  MICOM,  suggested 
that  a  very  effective  optical  system  for  missile  guidance  applications  would  be 
one  made  of  thin  films  of  mylar  or  similar  materials.  Such  a  mirror  or  system  of 
mirrors  provide  a  lightweight,  inexpensive  optical  element  or  system  of  elements 
and  it  has  the  potential  for  simple  refocussing.  By  simply  changing  the  pressure 
on  the  system,  one  could  change  the  curvature  of  the  mirror.  Other  possibilities 
for  shaping  the  mirror  exist.  For  example,  the  membrane  could  be  made 
aspherical  through  local  alteration  of  its  chemical  or  physical  structure.  These 
techniques  seemed  to  hold  great  promise,  so  an  exploratory  program  was 
undertaken.  This  program  included  the  preliminary  formation  of  mirrors,  design 
considerations,  and  prescriptions  for  a  Fresnel  mirror  to  ameliorate  certain 
difficulties. 

The  first  stage  of  the  research  program  comprised  an  extensive  literature 
search,  because  it  seemed  likely  that  this  concept  would  have  been  investigated 
previously.  Such  was  the  case,  but  unfortunately,  no  articles  were  pertinent. 
Applications  were  for  radar  and  solar  energy,  neither  which  require  the  surface 
perfection  required  by  visible  or  infrared  optical  systems.  Therefore,  we  were  on 
our  own.  We  continued  communications  with  Dr.  Gamble  who  had  carried  out 
other  exploratory  investigations. 

Original  tests  showed  that  reflectorized  mylar  can  be  formed  in  reasonable 
concave  and  convex  figures.  Pressures  ranged  near  atmospheric  -  approximately 
10%  above  or  below  one  atmosphere  would  produce  curvatures  that  can  be 
described  as  about  F/5.  Tests  were  made  with  a  2-in  diameter  ring  of  mylar. 
Higher  speeds  and  lower  F/numbers  introduced  "wrinklingn  at  the  mirror  edges. 
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This  is  not  unexpected,  because  the  surface  is  being  distorted  from  a  plane  to  an 
approximate  sphere.  The  wrinkling  is  a  function  of  membrane  elasticity.  It 
appears  that  to  obtain  very  fast  optics  will  require  surfaces  of  the  approximate 
final  shape  or  an  equivalent  technique. 

From  results  of  the  first  set  of  experiments,  and  through  further  discussions 
with  Dr.  Gamble,  we  designed  a  Fresnel  version  of  such  a  system.  The  Fresnel 
mirror  is  of  the  same  configuration  as  a  Fresnel  lens,  a  stepped  approximation  to 
the  curvature  desired.  The  system  was  designed  so  that  imaging  experiments  can 
be  carried  out  with  a  helium-neon  laser.  The  steps,  therefore,  must  be  in  half¬ 
multiples  of  0.6328  nm.  The  membrane  can  be  no  thicker  than  0.005  in  (- 125  nm). 
The  diameter  is  assumed  to  be  2  in  (-  50  mm).  An  F/3  system  would  have  a 
focal  length  of  150  mm  and,  therefore,  a  radius  of  curvature  of  300  mm.  The  sag 
then  is  1.04  mm.  Such  a  system  can  be  realized  in  1643  one-wavelength  steps  or 
larger.  The  largest  step  should  probably  be  about  one/tenth  the  thickness  of  the 
material,  or  about  12.5  pm,  giving  ten  steps  across  the  material.  A  reasonable 
compromise  is  20  radial  steps,  each  of  about  0.1  mm.  In  fact,  each  step  will  be 
158  wavelengths,  or  99.97  jim,  a  reasonable  step  height.  If  the  height  is  too 
shallow,  the  unavoidable  rounding  of  the  comers  hdS  a  deleterious  effect  on  the 
operation  of  the  Fresnel  element. 


PLANS 


The  Fresnel  master  has  been  ordered  from  the  University  of  Alabama, 
Huntsville,  and  will  be  delivered  to  Dr.  Gamble. 
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APPENDIX  A 

RADIOMETRIC  CALCULATIONS  AND 
MEASUREMENTS 


SOLAR  CONSTANT  -  JUSTIFICATION  OF  METHODS 


1.  Radiant  Exitance  of  sun  in  visible  spectrum,  from  Planck's  equation,  assuming  it 
to  be  a  perfect  black  body  at  5900  K: 


M  =  2.95  •  103  Watts/cm2 

2.  Radiance  of  sun,  assuming  that  it  is  a  Lambertian  source: 

L  -  —  L  *  940  Watts/cmVsr 

it 

3.  Solid  angle  subtended  by  sun  as  seen  from  earth: 

radius  of  sun  r  =  6.95  •  1010  cm 

distance  earth-to-sun  d  =  1.5  •  1013  cm 


a  -  6.7  •  10'sjr 


4.  Irradiance  incident  on  earth's  atmosphere: 

E  =  LQ  E  =  0.063  Watts/cm2 

E  •  104  =  634  Watts/m2 


5. 


For  comparison,  assuming  that  46%  of  solar  energy  falls  in  the  visible  band,  the 
solar  constant  is: 

0.46  •  1360  =  626  Watts/m2 


Conclusion:  The  irradiance  incident  on  earth's  atmosphere  as  derived  from  our 

radiometric  model  is  virtually  identical  to  the  value  given  for  die 
solar  constant  in  the  Manual  of  Remote  Sensing. 

Note:  353%  of  solar  energy  fells  in  .38  to  .66  nm  band. 


ADJUSTED  SOLAR  CONSTANT  AND  FILTER  TRANSMISSION 

Planck  Equation  for  Black  Body  Radiant  Exitance  of  sun: 

14388  1  -jT  (in  Watts/ cm2/ ^m) 

•  5900 J  J 

Curve  Fit  for  CS  1-56  infrared  blocking  filter: 

AO  :=  -23807.4  A1  :=  276308 

A2  :=  -1335780  A3  :=  3436050 

A4  :=  -4947200  A5  :=  3774250 

A6  :=  -1191400 


/  (X)  -  0.01  •  [40  +  AlX  +  A2X1  +  A3X*  +  A4X*  +  A5XS  +  A6X6]  (%  transmission) 


M(X)  •  37418  •  X-‘ 


5  [ex^ 


1. 


Calculated  Solar  Constant  at  sea  level,  adjusted  for  transmission  band  of  filter 

.66 


El 


0.83  •  6.7  •  10-5 


1 

~K 


•  J  M(k)  dX 

3t 


El  -  0.043  Watts/cm2 


2.  Average  filter  transmission  (experimental)  given  by  measured  value  divided  by 
Solar  Constant: 

77  -  -100  77  -  25.7% 

El 


3.  Calculated  Solar  Constant  at  sea  level,  adjusted  for  spectral  transmission  of  filter 

M 

E2  -  0.83  •  6.7  •  10-4  •  JL  •  j[/(X)  •  M(X)  dX 
E2  -  0.012  Watts! cm2 

4.  Average  filter  transmission  (calculated)  given  by  ratio  of  Solar  Constants  from 
parts  1  and  3: 


T2 


±1  -  100  72  -  28.5% 

El 
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WATTS/CM2  TO  LUX  CONVERSION  FACTOR 
Planck  Equation  lor  Black  Body  Radiant  Exitance  of  sun: 

M(k)  -  37415  •  X-5jexpj^il!E_J  -1 J*  (in  Watts/cm2/jim) 

Curve  Fit  for  photopic  response  curve: 

v(X)  «  1.0185  exp[-285.53  (k  -  .56017)*]  <rel-  response) 

Radiant  Exitance  of  sun  over  filter  transmission  band: 

M 

M-f  M(k)  dk  (in  Watts/cm2) 

.31 

Luminous  Exitance  of  sun  over  filter  transmission  band: 

M 

Mw  -  683  •  J  v(k)  M(k)  dk  (in  Lumen/cm2) 

31 

Conversion  Factor  (Watts  to  Lumen  or  Watts/m2  to  Lux): 

F  »  —  »  262  Lumen/Watt  or  Lux/Watt/m2 

K 

10*  •  —L  -  2.62  •  106  Lux/Watt/cm2 
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b  Transmission 


CONVERSION  FACTOR:  CD/M2  TO  PHOTONS/SEC-CM2-SR 


(sunlight  only  —  450  run  to  950  nm  band) 

Planck  Equation  for  Black  Body  Photon-Flux  Exitance  of  sun: 

Mp(k)  =  1.885  •  1023  •  X^Jexpj^JlHLj  -lj*  (in  photons/sec-onVpm 

Planck  Equation  for  Black  Body  Radiant  Exitance  of  sun: 

Mt(k)  »  37415  •  \'s  jexpj^ili^j  -lj  (in  Watts/cm^tm) 

Curve  Fit  for  photopic  response  curve: 

v(X)  -  1.0185  exp  [-285.53  •  (X  -  .56017)2]  (rel*  response) 

Photon-Flux  Radiance  of  sun  over  450  -  950  nm  band: 

Lp  -  iff  Af,(X)<*xl  (hi  photons /sec-cm2-sr) 


Luminance  of  sun: 


L  -  1  •  683  •  f  v(X)  MO.)  dk  (in  Cd/cm2) 

'it  J. 


Conversion  Factor  (Cd/m2  to  Photons / sec-cm2-sr) : 

F  -  •  lO*4  -  1.83  •  10“  Photons/sec-cm2-sr/Cd/m2 


CONVERSION  FACTOR:  CD/M*  TO  PHOTONS/SEC-CM*-SR 
(sunlight  only  -  450  nm  to  800  nm  band) 

Planck  Equation  for  Black  Body  Photon-Flux  Exitance  of  sun: 


M  (X)  =  1.885  ■  1023  • 

P 


exp  _  14388  -l  (in  photons/ sec-cm2/|im) 

X  •  5900 J  J 


Planck  Equation  for  Black  Body  Radiant  Exitance  of  sun: 

M,(K)  *  37415  •  X‘s  jexpj^J4^Lj  -lj  (in  Watts/cm2/pm) 

Curve  Fit  for  photopic  response  curve: 

v(X)  -  1.0185  exp  [-285.53  •  (X  -  .56017)2]  CreL  response) 

Photon-Flux  Radiance  of  sun  over  450  -  800  nm  band: 


Lp  «- 

p  it 


SO 

J  Mp(k)dk 

is 


(in  photons /sec-cm2-sr) 


Luminance  of  sun: 


.80 

L,  -  1  •  683  •  f  v(X)  Af,(X)  dX  (in  Cd/cm2) 
11  is 


Conversion  Factor  (Cd/m2  to  Photons/sec-cm2-sr): 


F  «  •  10“*  «  1.33  •  10“  Photons/sec-cm2-sr/Cd/m2 
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CONVERSION  FACTOR;  CD/M4  TO  PHOTONS/SEC-CM2-SR 
(sunlight  only  -  450  nm  to  670  run  band) 

Planck  Equation  for  Black  Body  Photon-Flux  Exitance  of  sun: 

Mp(X)  «  1.885  •  10“  •  XM  [^8»J  _ij  (in  photons/sec/cmVpm 


Planck  Equation  for  Black  Body  Radiant  Exitance  of  sun: 

M.(X)  =  37415  •  X's  |exp|^Jl|^_j  -lj  (in  Watts/cm2/pm) 

Curve  Fit  for  photopic  response  curve: 

v(X)  -  1.0185  exp  [-285.53  •  (X  -  .56017)2]  <*d.  response) 

Photon-Flux  Radiance  of  sun  over  450  -  670  nm  band: 


L'"7 


,S7 

f  Mp(X)dX 

is 


(in  photons /sec-cm2-sr) 


Luminance  of  sun: 


.67 

L,  -  1  •  683  •  f  v(X)  M'(X )  dk  (in  Cd/cm2) 
*  is 


Conversion  Factor  (Cd/m2  to  Photons/sec-cm2-sr): 


F  *  JL  •  10"*  -  8.41  •  10“  Photons/ sec-cm2-sr/ Cd/ m2 
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CONVERSION  FACTOR:  CD/M2  TO  WATTS/CM2-SR 
(sunlight  only  -  450  nm  to  950  nm  band) 

Planck  Equation  for  Black  Body  Radiant  Exitance  of  sun: 

M(K )  *  37415  •  X'5  -lj'  (in  Watts/cm2/nm) 

Curve  Fit  for  photopic  response  curve: 

v(X)  *  1.0185  exp  [-285.53  •  (X  -  .56017)2]  (rel.  response) 

Radiance  of  sun  over  450  -  950  nm  band: 

F95 

L  -_L  f*Af(X)  dk  (in  Watts/cm2-sr) 

*  Li 


Luminance  of  sun: 


.93 

L  -  1  •  683  •  f  v(X)  M(X)  dk  (in  Cd/cm2) 
jc  1. 


Conversion  Factor  (Cd/m2  to  Watts/cm2-sr): 

F  •  JL  •  10“*  -  5.47  •  10'7  Watts/ cm2-sr/Cd/m2 
L 


19 


CONVERSION  FACTOR:  CD/M2  TO  WAITS/CM2-SR 
(sunlight  only  -  450  nm  to  800  nm  band) 

Planck  Equation  for  Black  Body  Radiant  Exitance  of  sun: 

M(k)  -  37415  •  X'5  jexjjJjilHLj  -lj  (in  Watts/cm2/^m) 


Curve  Fit  for  photopic  response  curve: 

v(X)  =  1.0185  •  exp  [-285.53  •  (X  -  .56017)2]  (rel.  response) 

Radiance  of  sun  over  450  -  800  nm  band: 


W 


J*Af(X)  dk 
is 


(in  Watts/cm2-sr) 


Luminance  of  sim: 

to 

L.  -  1  •  683  •  f  v(X)  •  Af(X)  dk  (in  Cd/cm2) 

It  A€ 


Conversion  Factor  (Cd/m2  to  Watts/cm2-sr): 

F  •  h. 1  •  10~*  »  4.33  •  10*7  Watts/cm2-sr/Cd/m2 
L 
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CONVERSION  FACTOR:  CD/M2  TO  WAITS/CM2-SR 
(sunlight  only  -  450  nm  to  670  nm  band) 

Planck  Equation’fdr  Black  Body  Radiant  Exitance  of  sun: 

M(k)  -  37415  •  X*5  jexpj^llHLj  -lj‘  (in  Watts/cmVum) 


Curve  Rt  for  photopic  response  curve: 

v(X)  *  1.0185  exp  [-285.53  •  (X  -  .56017)2]  (rel.  response) 

Radiance  of  sun  over  450  -  670  nm  band: 

.St 

Lt  ■  _L  •  J  Af(X)  dk  (in  Watts /cm2-sr) 

Luminance  of  sun: 

St 

L  -  1  •  683  •  f  v(X)M(X)dX  (in  Cd/cm2) 
k  L 


Conversion  Factor  (Cd/m2  to  Watts/cm2-sr): 

L 

F  -  _i  •  10-*  -  3  •  10'7  Watts/cm2-sr/Cd/m2 
L 
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ESTIMATED  SIGNAL  FOR  MEASURED  TARGET  RADIANCES 


Object 
(in  full 
sunlight) 

Measured 
_  Value 
(E) 

(mW/cm1) 

Actual 

Value 

(E') 

(mW/cm1) 

Target 
Radiance  (L) 
(mW/cmVsr) 

Irrad.  on 
FPA 
(E*) 

(mW/cm2) 

Phot.  Flux 
Density 
(Ep) 

(photons/ 

sec-cm2) 

Signal 

(S) 

(electrons) 

Black  Roof 

0.065 

0232 

1.085 

0.192 

5.02E+14 

1.25E+Q5 

Red  Brick 

0.100 

0357 

1.669 

0295 

7.72E+14 

1.93E+05 

Dull  Aluminum 

0.150 

0536 

2503 

0.443 

1.16E+15 

2.89E+05 

Concrete 

0270 

0.964 

4506 

0.798 

2.08E+15 

5.21  E+05 

White  Roof 

0.290 

1.036 

4.840 

0.857 

224E+15 

559E+05 

Shiny 

Aluminum 

0.600 

2.143 

10.013 

1.772 

4.63E+15 

1.16E+06 

E  =  Irradiance  on  radiometer  head  (target  filling  entire  field  of  view) 

E'  =  E/0.28  =  Irradiance  on  radiometer  head  adjusted  for  average  filter  transmission 

L  =  E'/0.214  =  Target  Radiance  (0.214  =  solid  angle  subtended  by  target  =  field  of  view 
of  radiometer  head) 

Ed  =  (tc/4)  •  0.90AF/#)2  •  L  =  calculated  Irradiance  on  Focal  Plane,  assuming  F/2  lens 
with  90%  transmission 

Ep  =  2.612  •  1015  •  Ed  =  equivalent  Photon  Flux  Density  incident  on  Focal  Plane,  adjusted 
for  transmission  band  of  filter 


S  =  t\  •  t  -  A  •  Ep  =  number  of  electrons  generated  in  well,  assuming  a  Quantum 
Efficiency  On)  of  25%,  integration  time  (t)  of  1  ms  and  active  pixel  area  (A)  of  (10  pm)2 


MEASURED  VS.  EXPECTED  IRRADIANCE  VALUES 


Material 
(in  full  sunlight) 

Reflectance 

(P) 

Measured 
Irradiance  (Em) 
(mW/cm2) 

Expected 
Irradiance  (EJ 
(mW/cm2) 

Black  roof 

10% 

0.065 

0.075 

Concrete 

35% 

0.270 

0.262 

Aluminum 

90% 

0.600 

0.674 

E.  =  E,/jt  •  p  •  Q 


where 


E,  =  measured  solar  constant  (11  mW/cm2) 

p  =  reflectance  of  material  (rough  estimate  based  on  graphs  in  ER  Handbook) 
£2  =  solid  angle  of  field-of-view  of  radiometer  (0.214  sr) 
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APPENDIX  B 

EXPECTED  SIGNAL-TO-NOISE  RATIO 
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EXPECTED  SIGNAL-TO-NOISE  RATIO 

-  Signal  values  are  expected  to  vary  between  1E3  and  1E6  electrons  (see  Appendix  A). 

-  Noise  sources  are  photon,  read  and  dark  current  noise. 


-  Photon  noise:  np(S)  «  Js 

-  Read  noise:  nr  =*  100 

-  Dark  current  noise: 

active  area  A  =  225  •  10*  cm2 

integration  time  t  =  0.001  sec 

q  =  1.6  •  1  (T19  Coul. 

dark  current  Jd:  10*  A/cm2  at  300  °K 


/ 


nd  = 


Jd  •  A  •  t 


N 


<7 


nd  =  3.75 


-  Signal-to-Noise  ratio: 

SN  ( S,nr )  -  - S  __ 

Jnp(s)a+nr*+nd* 

-  S/N  vs.  signal  (log/log): 


1000 


SNlS.nr) 


9.5 


i  1  . 

V - « - 4 - -4 - -f  •  —  -4 - -4 -  -4  -  - 4 - 4  — 4 - * 

»'A !  1-f  -  U.<4«  .  *  '  ¥  -  *  «  -  •  i  *  %  « 

lOOO  S  1000000 


S/N  vs.  read  noise  (log/log): 

-  S/N  vs.  read  noise  (log/Iog): 


300 


.SEJtgJ&M-yJgy.. . .. _ ?  4-4W- 


SMS, nr)  } 


30  ^  ■  -  +; -  - „»  •  «  -  «  -« 

10  nr 


A  A  A  ■ 


1000 
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APPENDIX  C 
POWER  DISSIPATION 
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Power  Dissipation 


The  power  dissipation  on  a  focal  plane  array  depends  on  three  different  processes. 
Most  of  the  power  is  generated  by  capacitive  charge  and  discharge  currents  running  through 
resistive  elements  in  the  device  structure.  Another  significant  power  source  can  be  the  on- 
chip  output  amplifiers.  The  third  source  of  power  generation  are  carrier  lift  and  friction 
mechanisms  as  charge  packets  are  moved  across  the  device.  All  power  sources  depend  on 
basic  material  parameters,  the  device  architecture  and  layout,  and  the  clocking  scheme  used. 
The  follwoing  analysis  neglects  clock  pulse  degradation  due  to  lossy  clock  lines. 


As  a  clock  pulse  is  applied  to  the  device,  the  oxide  and  depletion  region  capacitances 
in  each  pixel  are  charged  up.  The  charge  current  dissipates  energy  by  flowing  through  the 
high-resistance  epitaxial  layer  into  the  substrate.  For  a  trapezoidal  clock  pulse  with  rise  and 
fall  times  tr  and  voltage  differential  Vf  between  the  two  clock  bias  levels,  the  current  is 
approximately  constant  over  the  linear  voltage  ramp  and  can  be  expressed  as 


C_.K 


The  dissipated  power  is  then  given  by  integrating  PR  over  the  period  of  one  clock  cycle. 
Since  the  current  is  zero  except  when  the  clock  voltage  is  changing,  this  integral  reduces  to 


where  T  is  the  period  of  the  clock  cycle  (T=  l//c).  The  factor  of  two  is  due  to  the  power 
being  dissipated  on  both  the  rising  and  the  falling  edge  of  the  clock  pulse.  To  get  the  total 
power,  P,pi  has  to  be  multiplied  by  the  number  of  affected  pixels  and  the  duty  cycle  of  the 
clock  waveform.  Solving  the  integral, 


=  2  , 

where  #  is  the  number  of  phases,  N  is  the  number  of  pixels  and  d  is  the  duty  cycle. 


A  second  source  of  power  dissipation  is  the  resistance  of  the  polysilicon  lines.  This 
resistance  is  encountered  as  the  the  current  charges  up  the  pixel  capacitance  Cpix  and  the 
poly-to-poly  capacitance  C^.  Following  the  above  argument, 


This  current  has  to  be  summed  over  half  of  each  1024-pixel  clock  line,  since  each  pixel  sees 
a  different  current.  The  total  power  dissipated  in  the  polysilicon  clock  lines  is 


A-l 
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Here  N{  is  the  number  of  lines  and  zn2=m(m+  l)(2m+ 1)/6.  The  additional  factor  of  two 
is  due  to  the  fact  that  each  line  is  1024  pixels  long  and  driven  on  both  ends,  so  that  m*=512. 


Power  is.  also  generated  due  to  carrier  lift  and  friction  as  charge  packets  are 
transported  across  the  array.  This  power  is  given  by 

*m  -  N~mld 


S  is  the  number  of  electrons  in  each  signal  package,  L  is  the  pixel  length,  q  the  electron 
charge  and  n  the  carrier  mobility. 


Finally,  power  is  generated  by  the  64  output  FETs  on  the  chip.  Neglecting  rise  and 
fall  times  on  the  output  signal,  it  is  given  by 

=  «(»%>- 

VOD  is  the  output  drain  supply  voltage,  Va  is  the  output  voltage  and  Rlomd  is  the  load 
resistance  attached  to  the  output  amplifier. 


These  equations  can  now  be  evaluated  for  a  given  array.  In  this  case  the  Ford 
F1024B  array  running  in  frame  transfer  mode  at  a  frame  rate  of  1000  fps  will  be  used.  To 
calculate  the  total  power,  the  readout  process  has  to  be  separated  into  three  sections.  First, 
the  charge  is  transferred  from  the  integration  region  to  the  storage  region  at  a  clock  speed 
of  1  MHz.  This  process  takes  256  ns.  During  the  remaining  744  /xs  of  the  frame  period,  the 
256  lines  of  data  in  the  storage  region  are  read  out  into  the  serial  registers  at  a  clock  speed 
of  344  kHz.  At  the  same  time,  the  serial  registers  are  clocking  out  the  charge  at  16.8  MHz, 
running  idle  for  1  ns  every  32  cycles  to  wait  for  a  new  line  of  data  to  be  transferred  from  the 
storage  area. 


A  few  special  rules  apply.  The  outputs  of  the  device  are  at  a  low  signal  level, 
corresponding  to  the  full  well  output,  for  half  the  time  during  which  the  serial  registers  are 
running  (0.244  duly  cyle).  For  the  remaining  time,  they  are  at  the  reset  level  V0\  which  is 
closer  to  For  the  serial  readout,  the  pixel  capacitance  is  doubled  and  the  epi  resistance 
halved  since  the  pixels  in  the  serial  readout  register  have  twice  the  area  of  the  array  pixels. 
Finally,  the  number  of  pixels  for  the  active-to-storage  area  transfer  for  the  power  generated 
by  carrier  lift  and  friction  is  halved  since  only  half  the  pixels  on  the  array  carry  charge. 
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Epitaxial  Laver 


Transfer  to  Transfer  to  Serial  Register 

Storage  Area  Serial  Register  Readout 


4. 

3 

3 

' 

3 

N 

1048576 

1048576 

2048 

Q»ix 

1.5 

fF 

1.5 

fF 

3 

fF 

R.pi 

80 

kn 

80 

kn 

40 

kn 

Vf 

10 

V 

10 

V 

12 

V 

t. 

100 

ns 

100 

ns 

10 

ns 

U 

1 

MHz 

344 

kHz 

16.8 

MHz 

d 

25.6 

% 

74.4 

% 

48.8 

% 

P.pi 

290 

nWatts 

150 

nWatts 

520 

nWatts 

total: 

960 

nWatts 

Polvsilicon  Lines: 

Transfer  to 
Storage  Area 

Transfer  to 

Serial  Register 

Serial  Register 
Readout 

* 

3 

3 

3 

N, 

1024 

512 

2048 

Sn* 

44870400 

44870400 

1 

&ix 

1.5 

fF 

1.5 

fF 

3 

fF 

Cpoly 

10 

fF 

10 

fF 

20 

fF 

Rpo  ,y 

100 

n 

100 

n 

200 

n 

V, 

10 

V 

10 

V 

12 

V 

tr 

100 

ns 

100 

ns 

10 

ns 

U 

1 

MHz 

344 

kHz 

16.8 

MHz 

d 

25.6 

% 

74.4 

% 

48.8 

% 

P.p. 

6.5 

Watts 

3.3 

Watts 

540 

/i  Watts 

total: 

9.8 

Watts 

Output  fr'ETs: 

_ YJgw  V,  high 


Voo 

20 

V 

20 

V 

11 

V 

12 

V 

10 

kn 

10 

kn 

d 

24.4 

% 

75.6 

% 

PfET 

160 

mWatts 

460 

mWatts 

total:  620  mWatts 
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Carrier  Lift/Friction: 


Transfer  to 

Transfer  to 

Serial  Register 

Storage  Area 

Serial  Register 

Readout 

*_ 

3 

3 

3 

N 

524288 

262144 

1024 

Cox 

26 

fF 

26 

fF 

52 

fF 

S 

100000 

e* 

100000 

e' 

100000 

e* 

L 

15 

fim 

15 

/itn 

15 

Mm 

M 

1300 

cm2 /Vs 

1300 

cm2 /Vs 

1300 

cm2 /Vs 

v. 

10 

V 

10 

V 

12 

V 

fc 

1 

MHz 

344 

kHz 

16.8 

MHz 

d 

25.6 

% 

74.4 

% 

48.8 

% 

art 

62 

mWatts 

31 

mWatts 

4.8 

mWatts 

Pfric 

3.7 

MWatts 

0.6 

MWatts 

3.9 

MWatts 

total: 

98 

mWatts 

grand  total: 

10.5 

Watts 

The  total  power  dissipated  by  the  array  in  this  model  is  10.5  Watts.  This  points  out 
several  problems.  The  main  power  contribution  of  9.8  Watts  is  generated  in  the  array 
polysilicon  lines  due  to  the  extremely  high  capacitances  encountered  between  adjacent  lines. 
In  an  actual  device,  the  power  dissipation  would  never  reach  these  levels  since  the 
assumption  of  lossless  clock  lines  does  not  hold.  The  clock  pulse  would  be  degraded 
significantly  towards  the  center  of  a  clock  line,  so  that  less  current  is  generated  and  the 
power  dissipation  is  kept  low. 

Two  solutions  exist  for  this  problem.  The  polysilicon  clock  lines  could  be  strapped 
by  aluminum  lines,  which  have  less  resistance  and  therefore  cause  less  pulse  degradation  and 
power  generation.  However,  since  the  large  capacitances  remain,  this  solution  would  create 
extremely  high  currents  during  the  rising  and  falling  edges  of  the  clock  pulses.  Measures 
have  to  be  taken  to  reduce  the  line-to-line  capacitance,  either  by  reducing  the  overlap  area 
or  by  increasing  the  thickness  of  the  oxide  insulation  in  between  polysilicon  layers. 
Calculations  show  that  by  decreasing  the  overlap  area  by  half  to  one  micron  overlap,  and 
increasing  the  aside  thickness  from  1000  A  to  5000  A,  the  overlap  capacitance  can  be 
reduced  to  1  fF.  This  in  turn  reduces  the  total  dissipated  power  to  one  Watt  The  best 
solution  is  probably  a  combination  of  both  aluminum  lines  to  ensure  accurate  clock  pulses 
and  reduced  capacitance  to  reduce  high  charging  currents. 
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APPENDIX  D 

CONTRAST  DEGRADATION  DUE  TO 
INTEGRATION/READOUT  TIME  RATIO 
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Contrast  Degradation  due  to  Integration/Readout  Time  Ratio 

(T.  Graeve) 

Let  the  ratio  of  integration  time  to  readout  time  be  R  (i.e.  R:l).  For  a  given  array 
size  and  frame  rate,  a  full  frame  architecture  has  the  lowest  R,  a  frame  transfer  device 
has  a  higher  R,  and  an  interline  device  has  the  highest  R.  However,  for  the  latter  there 
is  no  degradation  in  contrast  due  to  charge  readout.  The  following  analysis  is  therefore 
valid  only  for  full  frame  and  frame  transfer  devices. 

For  a  normal  scene,  we  define  the  maximum  and  minimum  irradiance  incident 
on  the  devices  EmiX  and  E^  respectively.  The  maximum  contrast  for  this  scene  then 
becomes 

f*  _  ^inm) 

+  2 t,n  • 


where  rj  is  the  quantum  efficiency  and  the  area  of  the  device.  During  readout,  each 
pixel  collects  an  additional  charge  proportional  to  the  average  irradiance  on  the  column 
that  it  is  transferred  through.  If  this  average  irradiance  is  given  by  a  •  E,^  where  a 
is  a  factor  between  zero  and  one  and  the  readout  time  is  t,/R,  the  minimum  and 
maximum  signal  levels  are 


in  general,  the  reduced  formula  can  be  applied  as  long  as  the  dynamic  range  of  the 
incurrent  radiation  is  greater  than  20  dB  and  the  dark  current  in  the  device  is  less  than 
100  nA/cm2. 
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Note  that  this  analysis  assumes  that  all  random  noise  sources  have  been  averaged 
to  zero.  The  remaining  noise  sources  that  average  to  a  DC  level  are  included  in  the  dark 
current  term.  The  resulting  contrast  degradation  is  superimpressed  onto  the  MTF 
response  of  the  device  (i.e.,  it  scales  down  the  MTF  curve). 

Example:  Under  the  following  conditions  (either  equation  applies)  the  contrast  becomes: 

Ew  =  5  •  1014  photons /sec- cm2  q  =  25% 

Emin  =  1  •  1012  photons /sec-cm2  Ad  =  (15  nm)2 

n  =  JdAd/q  =  1.4  •  104  electrons/sec 

(for  a  dark  current  of  InA/cm2) 


a  =  1  and  R  =  17 

=> 

C  =  89% 

a  =  1  and  R  =  3 

=> 

C  =  60% 

a  =  1/2  and  R  =  17 

C  =  94% 

a  =  1/2  and  R  =  3 

=> 

C  =  75% 
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contrast 
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APPENDIX  E 
SPICE  SIMULATIONS 


Analysis  from  1/30/91) 


Volts  Amps  (IE-4) 


APPENDIX  F 

MULTISPECTRAL  RADIOMETER 
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narrow  FOV 


wide  FOV 


Multispectral  Radiometer 

dimensions: 

weight: 

power  requirements; 
operating  temp.: 
video  output: 
data  output 

camera: 

WFOV  radiometer 

NFOV  radiometer: 


data  capture  rate: 
spectral  Biters: 


System  Specifications 

24"  (L)  x  8"  (W)  x  12"  (H) 

80  lbs. 

28  Vdc,  2  amps 
-15  °C  to  40  °C 

EIA  standard,  see  Xybion  camera  specs 

8-bit  parallel  interface  to  Datum  Video 
Data  Encoder 

30.8°  (H)  x  23.3°  (V)  field  of  view 
16mm,  F/1.4  lens 
24.2°  circular  field  of  view 
50  mm,  F/1.8  lens 

dynamic  range  from  1  nW/cm2-sr  to  0.85 
mW/cm2-sr  extended  source  radiance  at 
maximum  iris,  depending  on  filter 

average  standard  deviation  on  50 
samples  is  9.6  D.N. 

4°  circular  field  of  view 
135  mm,  F/2.8  lens 

dynamic  range  from  50  nW/cm2-sr  to  2 
mW/cm2-sr  extended  source  radiance  at 
maximum  iris,  depending  on  filter 
average  standard  deviation  on  50 
samples  is  11.5  D.N. 


60  samples/sec 


filter  # 

wavelength 

bandwidth 

1 

400  run 

70  nm 

2 

450  nm 

80  nm 

3 

500  nm 

80  nm 

4 

600  nm 

80  nm 

5 

700  nm 

80  nm 

53 


